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Abstract 

The quark form factor is known to exponentiate within the framework 
of dimensionally regularized perturbative QCD. The logarithm of the form 
factor is expressed in terms of integrals over the scale of the running coupling. 
I show that these integrals can be evaluated explicitly and expressed in terms 
of renormalization group invariant analytic functions of the coupling and of 
the space-time dimension, to any order in renormalized perturbation theory. 
Explicit expressions are given up two loops. To this order, all the infrared and 
collinear singularities in the logarithm of the form factor resum to a single 
pole in e, whose residue is determined at one loop, plus powers of logarithms 
of e. This behavior is conjectured to extend to all loops. 
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1 Introduction 



The asymptotic behavior of the electromagnetic form factor of charged particles at 
high energy has been the object of theoretical studies for almost half a century. It 
is the simplest gauge theory amplitude to exhibit a doubly logarithmic behavior, as 
a consequence of its infrared and collinear singularities. Thus the understanding of 
its asymptotic energy dependence requires a resummation of perturbation theory, 
going beyond renormalization group logarithms. This resummation was performed 
for the first time, at the leading logarithm level, by Sudakov [|TJ, who condidered 
the off-shell form factor for an abelian gauge theory. He found that the leading 
(double) logarithms exponentiate, which results in the strong suppression of the 
elastic scattering of charged particles at high energy which bears his name. A 
similar exponentiation for the on-shell form factor was obtained in , and in later 
years the result was extended to nonabelian gauge groups 0. 

Further extending the exponentiation to all subleading logarithms is a nontrivial 
task, which requires the complex machinery of perturbative factorization theorems. 
This goal was achieved by Mueller [[| and Collins || for abelian theories, and shortly 
later by Sen |J for QCD. In essence, the complete exponentiation is achieved by 
showing that the energy dependence of the form factor must follow an evolution 
equation, which embodies the constraints of renormalization group and gauge in- 
variance. This evolution equation was typically solved in terms of an undetermined 
initial condition, containing low-energy information and depending on some chosen 
infrared regulator, say a gluon mass. It was later realized [J] that dimensional regu- 
larization can be used directly at the level of the evolution equation. One can then 
give a simple and explicit solution for the form factor, using as initial condition 
the electric charge at vanishing momentum transfer. When the solution is written 
in this form, all infrared and collinear poles are explicitly exponentiated, and the 
results can be directly compared with diagrammatic calculations. 

The quark form factor, although divergent, is a quantity of considerable in- 
terest for phenomenology: in fact, it appears in the computation of several cross 
sections, and techniques similar to the ones used for its exponentiation have been 
applied to several processes of interest in high energy physics, in particular to pro- 
cesses in which the perturbative series needs to be resummed because of logarithmic 
enhancements in special regions of phase space. Outstanding examples are the in- 
clusive Drell-Yan cross section || and the transverse momentum distribution of 
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vector bosons produced in hadronic collisions |J. More recently, similar techniques 
have been generalized to hadronic cross section with colored particles in the final 



state ||10fl . Strikingly, the quark form factor enters directly in the hard partonic 
cross section for the Drell-Yan process in the DIS scheme, which is proportional to 
the modulus squared of the ratio of the timelike to the spacelike form factor [Rl O . 



In f?j it was shown that this ratio is given by an infinite phase (a generalization of 
the Coulomb phase, expressed by a series of pure counterterms in the MS scheme), 
times a manifestly finite exponential factor. 

In the present paper, I will build upon the results of Ref. J7|, and show that 
the resummed expression for the form factor can be explicitly evaluated to the 
desired accuracy, thus expressing the logarithm of the form factor as an analytic 
function of a s (Q 2 ) and e, manifestly renormalization group invariant. The key 
observation is the following: typically, in resummed expressions, integrals over the 
renormalization scale cannot be explicitly evaluated, because the one-loop running 
coupling has a Landau pole on the integration contour, which extends from the 
hard scale, /i 2 = Q 2 , all the way down to fi 2 = 0, in the region dominated by 
nonperturbative effects. The sensitivity to the Landau pole, which can be estimated 



using various methods and models |T2| , is then taken as a measure of the influence 
of nonperturbative effects on the given observable. In dimensional regularization, 
however, the running coupling depends not only on the renormalization scale, but 
also on the dimensional regularization parameter e. At one loop, the coupling still 
exhibits a Landau pole, but for general (complex) e the pole is located at complex 
values of the renormalization scale fi 2 , away from the integration contour. As a 
consequence, all integrals appearing in the logarithm of the form factor are well- 
defined, and can be evaluated. Alternatively, one can express the one-loop running 
coupling as a power series in terms of a s (Q 2 ), and integrate the series term by term: 
one finds that there is no factorial growth in the coefficients of the series obtained 
after integration; in fact, the series can be summed to reconstruct the same analytic 
function of a s and e that was found by direct integration. It is easy to see how these 
results generalize when the running coupling is expressed including higher orders in 
the QCD f3 function. The complexity of the answer grows, but the most relevant 
features remain unchanged; in particular, the answer is manifestly invariant under 
changes in the renormalization scale, to the given accuracy in the (3 function; also, 
infrared and collinear poles resum, so that the logarithm of the form factor has only 
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a single pole in e, as well as cut on the negative real axis in the e plane; conventional 
perturbative results can be recovered by reexpanding in powers of a s for finite e. 

The outline of the paper is the following: in Section |2|, I briefly review the known 
results on the exponentiation of the quark form factor. In Section ^ I illustrate 
the power of dimensional regularization by solving explicitly to all orders the re- 
cursion relation arising from the renormalization group equation for the generalized 
Coulomb phase: all the towers of poles (leading, next-to-leading, etc.) arising in 
the series of counterterms can be expressed in terms of the perturbative coefficients 
of the relevant anomalous dimension, and of the QCD (3 function; all the resulting 
series can be summed, and they give functions of e with at most a logarithmic sin- 
gularity; the leading behavior as e — > can be explicitly computed to all orders. In 
Section |], I turn to the main issue, and compute the logarithm of the form factor 
using one-loop results only, as well as the one-loop running coupling; the result is 
remarkably simple, a single dilogarithm and a single logarithm of functions of a s 
and e, explicitly independent of the choice of renormalization scale. In Section I 
generalize the calculation to two loops, obtaining again a relatively simple answer, 
containing only logarithms and dilogarithms. It is fairly clear that the calculation 
could be formally pushed to even higher orders, adding to the complexity of the 
answer, but without changing its basic features; furthermore, the perturbative co- 
efficients of the various functions appearing in the logarithm of the form factor are 
known only to two loops. Section |] summarizes the results and briefly discusses 
possible applications and future developments. 



2 Exponentiation of the form factor 

The techniques needed to compute the quark form factor in perturbative QCD to 
all orders in the logarithms of the energy are described in detail in Ref. Ul3| . Here 
I will just give the relevant definitions and sketch the main results, in the spirit of 
Ref. 0. 

The (timelike) quark form factor in dimensionally regulated massless QCD is 
defined by 

r M ( Pl ,p 2 ;/i 2 ,e) = (0\U0)\pu P 2) (2.1) 



Q 
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-iee q v(p 2 )%u(j) l ) T [—,a s (fi ),e 
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where J^{x) is the electromagnetic quark current, while \p\1P2) is a state containing 
a quark and an antiquark with momenta p\ and P2, and total energy Q 2 = (pi +P2) 2 , 
taken to be much larger than the confinement scale. In Eq. ( |2.1| ) I made use of 
the fact that in the massless limit the perturbative form factor has only one spin 
structure, and I defined the scalar form factor F to be normalized to 1 in the absence 
of strong interactions. Eq. ( |2.1|) is written in the renormalized theory, so that all 
ultraviolet divergences arising in perturbation theory have been dealt with, and one 
can take e = 2 — d/2 < to regulate infrared and collinear poles. Notice that, as a 
consequence of electromagnetic current conservation, the anomalous dimension of 
the form factor vanishes. 

The conventional analysis of the form factor Jl3| proceeds through the steps that 



are characteristic of the derivation of perturbative factorization theorems: one starts 
by identifying the integration regions in momentum space that contribute to the 
form factor at leading powerQ of Q 2 ; next, one constructs a factorized expression for 
the form factor, in terms of (nonlocal) operators, typically defined in terms of eikonal 
lines, whose matrix elements reproduce the leading behavior in the various regions. 
Each operator factor will typically depend both on the gauge choice and on the 
renormalization scale; imposing gauge and renormalization group (RG) invariance 
leads to an evolution equation describing the energy dependence of the whole form 
factor. In dimensional regularization, this equation takes the form 








Q2 

K (e, a s (> 2 )) + G — , a s {^ 



• (2-2) 

The functions K and G are characterized as follows: K is a pure counterterm, thus 



a series of poles in the MS scheme, which I will use throughout; G contains all the 
energy dependence, and is finite as e — > 0; the sum K + G must be renormalization 
group invariant, because the form factor is; thus K and G renormalize additively, 
according to 



( d d \ 

= - I H-q- +/3(e,a s )— 1 K(e,a s ) . 

technically, this is done using a massive regulator; in dimensional regularization all the energy 
dependence is logarithmic, leading regions are responsible for infrared and collinear poles, and one 
is in fact working with the whole form factor. 
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The anomalous dimension 7^ is well known in perturbative QCD. It is essentially 
the anomalous dimension associated with a cusp in a Wilson line, corresponding 
here to the hard vertex where the quark and the antiquark annihilate ||14|| . It is 
independent of e, and can be expanded perturbatively as 



7*(« s ) = E7l?R 



n=l 



a, 



71 



with coefficients known to two loops ]T5] 



7? 



2C, 



(2) 

Ik 



67 



-t(2))C A C F --n f C F 



(2.4) 



(2.5) 



Since we are dealing with divergent quantities, it is important to keep e < 
throughout. Thus the j3 function is defined by 



(3(a s 



dfx 



-2ea s + (3(a s ) 



a 



2 00 



Z,H n=0 



(2.6) 



At the leading nontrivial order, the solution of Eq. (|2.6f) is 



a ( ^2,a s (/Xo),e ) = 



^0, 



1 



A, 



Ait 



(2.7) 



with b = {11C a — 2n/)/3. Using this form of the running coupling, and its higher 
order generalizations, one can solve the renormalization group equation for the 
function G, as 



G 



^ 2 



G 



_( -Q 2 
1, a I — — , a s (//) 



+ r 



-Q 2 A 



dA 2 / /A 2 a 
7x [a — , a s (/x J 



(2.8) 



where the initial condition was chosen to emphasize that G is real for negative Q 2 . 

The second important consequence of dimensional regularization is the fact that 
it provides a natural initial condition for Eq. (|2.2j ), in terms of which one can 
construct an explicit solution. In fact, one readily recognizes that each term in the 
perturbative expansion of T(Q 2 ) must be proportional to a positive integer power 
of (fj, 2 1 \—Q 2 )) e . Thus for fixed e < all perturbative corrections to the form factor 
vanish in the limit Q 2 — > 0, and one can use 



r (0,a s (yU 2 



1 . 



(2.9) 
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One finds then 

„(Q 2 



+ G 



2 - 



a s (fi 



exp 



K (e, a s (/i 2 



-l,a { -^,a s (/i 2 ) 



1 r-Q 2 de 

2 Jo ~e 

1 r^rfA 2 / /A 2 , 2 



(2.10) 



Eq. (|2.10|) can be used directly as a generating function of diagrammatic results 
in dimensional regularization; for example, using one-loop results for the various 
functions appearing in the exponent, as well as the one-loop running coupling re- 
expanded in terms of a s (Q 2 ), one reproduces correctly the quartic and cubic poles in 
e at 0(a 2 ). Introducing the two-loop contribution to the anomalous dimension 7^, 
one matches the double pole as well; the only singularity requiring a full two-loop 
computation (the knowledge of G^ 2 \ as well as 0(e 2 ) terms in G^), is the single 
pole, in agreement with the general considerations of Ref. [fL6j . In Ref. |7| it was 
shown that Eq. ( |2.10| ) can be used to evaluate explicitly the ratio of the timelike 
to the spacelike form factor, in terms of a contour integral in the complex plane of 
the renormalization scale £ 2 . The result (at /i 2 = Q 2 ) is of the form 

T(l,a s (Q 2 ),e) 
r(-l,a s (Q 2 ),e) 



exp (e, a s (Q 2 )) + i £ d6 G (-1, a (e ie , a s (Q 2 ), e) , e 

X - / e rf0 7 x(«(e i0 ,a s (Q 2 ), e ))|} . (2.11) 



Eq. ( 2.11 ) proves that the ratio of form factors is given by an infinite phase times 
a finite exponential factor. The derivation of Eq. ( j2.11|) also shows that the ratio is 
dominated by ultraviolet contributions: in fact, it is expressed by a contour integral 
that never approaches the origin. Thus, one may conclude that it is free of power 
corrections related to the Landau pole. In particular, the integrals in Eq. ( 2.11|) can 
be evaluated explicitly to the desired accuracy in the running coupling; the result 
is of direct phenomenological relevance, since the modulus squared of the ratio of 
form factors appears in the resummed partonic Drell-Yan cross section 
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I will now concentrate on Eq. ( 2.10Q , and examine the structure of the exponent 
in more detail. In particular, I will be concerned with the fate of the Landau pole, 
that is expected to make the integration over the infrared region ill-defined. It will 
become apparent that dimensional regularization provides a solution to this prob- 
lem. I will however start with an exercise in the application of the renormalization 
group, which displays in a simple context the power of dimensional regularization. 
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3 The counterterm function 



Let us begin by considering the "Coulomb phase" K(e,a s ). In any minimal renor- 
malization scheme, it is a series of poles which can be written asf] 



00 K (a ) 

K(e,a s ) = £^f< (3.1) 

n=l 6 

00 , 

K n (a s ) = £^i m) (v 



In such a scheme, K(e, a s ) has no explicit dependence on the renormalization scale 
/i. The renormalization group equation (|2.3|) then reduces to 



d 

/3(e,a s )-—K(e,a s ) = -j K (oc B ) , (3.2) 

which, given the finiteness of 7^ and the form of /3(e,a s ), turns into a recursion 
relation for the perturbative coefficients of K. Explicitly 

as^—K^ots) = \^ K {a s ) , (3.3) 
aa s 2 

a s ^—K n+1 (a s ) = \(3{a s )—K n (a s ) . 
aa s l aa s 

Eq. ( p.3[ ) expresses the fact that the physics of the function K(e, a s ) is contained in 
the residue of the simple pole, Ki(a s ), which is completely determined by the 
anomalous dimension 7^, and in turn determines the coefficients of all higher 
poles. One can approximate the exact recursion relation by truncating the four- 
dimensional j3 function (3(a s ) to some chosen order. 

It is useful to collect the terms in K corresponding to leading poles ((a s /e) n ), 
next-to-leading poles (a" +1 /e n ), and so on. This is done by writing 



K(e,a s ) = £/C m (e,a s ), (3.4) 

m=0 

00 



m=0 

. / a \ n+m 1 



n=l 



7T / e 



According to general arguments, one expects the leading poles in K to be completely 
determined to all orders by the one-loop anomalous dimension and the one-loop (3 



2 Notice that K(e, a s ) has only one pole per loop. The leading pole in the exponent of T(Q 2 ), 
(a"/e" +1 ), is generated at each order by the integration of the anomalous dimension jk- 
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function; similarly, the m-th tower of poles, K, m is expected to require the anomalous 
dimension and the (3 function tom + 1 loops. 

In the remainder of this section, the following facts about these towers of poles 
will be established: 

• the recursion relation for the perturbative coefficients can be explicitly 
solved including all orders in the j3 function. 

• All the resulting series of poles, JC m (e, a s ), can be summed. 

• Upon summation, JC m (e,a s ) is an analytic function of a s and e, regular as 
e — > for m > 0. The only singularity at e — > is logarithmic and completely 
determined by a one loop calculation. 

• The finite limits /C m (0, a s ) (m > 0) can be computed for all m in terms of the 
perturbative coefficients of (3 and 7#. They reconstitute a power series in a s . 

As with most recursion relations, it is useful to start with the simplest cases. 



3.1 One loop 

The first of Eqs. ( fO| ) is readily solved and yields 



K ( r ] = ^iP ■ (3-5) 

In the second of Eqs. ( p.3| ) one can start by truncating the (3 function at 0(a%), 
keeping only the one-loop coefficient 60 • Using Eq. ( |3.5| ) one readily finds an ex- 
pression for the generic coefficient K^, 

which, as we will verify, is exact for n = m. Then 

Notice that this implies that the resummation of all the leading poles in the ratio 
( 12.111 ) degrades the singularity to a logarithm: 



r(l,a s (Q 2 ),e) 



r(-l,a,(QV) 



*>0 



1 + ^1 . (3, 



L.P. V 4?Fe 
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Since all other poles are perturbatively weaker, one might expect that they should 
resum to even weaker singularities, or to regular functions. This in fact turns out 
to be the case. 



3.2 Two loops 



To resum next-to-leading poles, one must include the next coefficient of the (3 
function, bi, as well as 7^'. Eq. ( |3~3|) may still be solved, and gives 



prim) _ _L 
n " 2m 



"-' 1 n-l 
k 



E 

k=0 



' 4 



ra— 1 — k 



bl \ (m-n-k+l) 



(3.9) 



where the sum is further truncated by the requirement that 7^ = for p < 1. 
One recovers the previous answer, Eq. ( |3.7|) , for the leading poles. Next-to-leading 
poles are given by 

n— 3 r , 



Kn - 1 ~2n[ 4 



bo 



^-(n-2)^ 



4 



(3.10) 



Again, the corresponding series is summable, and yields 



71=1 



71+1 I 



k / 27^6! 



12 



7T 



60 

2 (7^ 

&0 



7 {? I In I 1 + 



boa s 
Aire 



(3.11) 



27^61 



Aire 



bo b a s + 47re 



The logarithmic singularity receives a new contribution, which, however, is su- 
pressed by a power of e. Thus /Ci has a finite limit as e — > 0, given by 



£i(0, a s ) = 

7TOn 



(2) 

7x 7 



(3.12) 



3.3 All loops 



Encouraged by the simple form of the answers obtained for /C and JC%, one sets out 
to include higher orders in (3(a s ). It turns out that the recursion relation Eq. 
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is solvable to all orders, and the solution is a straightforward generalization of Eqs. 
~M and (B]l). One finds 



1 



2m ^ 



pi,...,p„_i=0 



n-1 



J'=l 



( m _ n+ l_£™ =1 p.) 

7x 



(3.13) 



To identify explicitly the towers of subleading poles, it is useful to reorganize the 
sum isolating the perturbative coefficients of the anomalous dimension jk. This is 
done by writing Eq. fl3.13Q as 



Qn-2 



1 CO *?i 

*i m) = i- E 7r n+1 ~ 9l) E ••• E 

Zm <?i=0 I 92=0 g„-i=0 



n-1 

n 

j=i 



4 



(3.14) 



where the sum over q\ is effectively cut off by the fact that for qi > m — n + 1 
the summand vanishes. It is not trivial to reconstruct the general form of the m-th 
tower of subleading poles from Eq. (|3.14j) ; however, upon generating a sufficient 
number of examples, one recognizes that 



/T (n) = — I — — 

- Ii n— m r, " 

In 



bo 
4 



n-m-l m+ i 
P =i 



J2C p (n,m)4 m+2 - p) 



(3.15) 



where 



C p (n, m) = E 

fc=0 



n — m — 1 



h, ■ ■ ■ , bp-i) 



(3.16) 



The coefficients C k (&i, . . . , &p-i) are polynomial expressions in the fej's, given by 



4 p) (&i, 



p-i 

E 

m» = 1 



(3.17) 



The key observation concerning Eqs. (|3.15|) - (|3.17|) is that the dependence on n is 
simple, so that summation over n can explicitly be performed for all m. In fact, for 
fixed m, Eq. (|3.15 ) can be rewritten as 



(n) 



1 

2n 



n 



(3.18) 



r=0 



with easily determined coefficients T> r (m). All the series to be summed are thus 
simple power series, in fact derivatives of the logarithmic series which resums the 
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leading poles, up to subtractions of a finite number of terms. It is easy to verify 
that the results obtained for fC m (e,a s ) are finite in the limit e — > 0. One can in 
fact determine the form of this limit explicitly for all m. It is given by polinomial 
expressions in the coefficients of the (3 function, similar to Eqs. ( |3.15|) - (|3.17|) : 



/C m (0, a s 



- ( — ) f E B p (m) - 1 



(m+l—p) 



where 



and 



j=0 



bo 



Bj{m,p) 



E 

q% = l 



n 

i=l 



(3.19) 



(3.20) 



(3.21) 



The most amusing fact here is perhaps that in the limit e — > one recovers a finite 
perturbative expansion in powers of a s /rt. 

The conclusion of this exercise is that the counterterm function K(e, a s ) is com- 
pletely determined, as an analytic function of e, a s , and the perturbative coefficients 
of j3 and jk- in the limit e —>■ this function has a logarithmic singularity, arising 
from the resummation of the leading poles, and a contribution which is indepen- 
dent of e and can be computed as a perturbative expansion in powers of a s /ir, with 
coefficients determined again by the functions j3 and jx- Explicitly 



K(e, a s ) = K DIV (e,a s ) + K FIN (a s ) + 0(e, elne) 



with K DIV (e, a s ) given by Eq. (|3.7|) , and 

oo 

K FIN (a s ) = ^/C m (0, a 



m=l 
in 



bo 



E 

m=l 



I /a 



m \ 7r 



.A 



(3.22) 



(3.23) 



p=0 



None of the series involved in this computation shows signs of factorial growth in 
the coefficients. In the final expression, Eq. (|3.23|) , the growth of the coefficients is 
implicit in their dependence on bi and 7^, the series in Eq. ( |3.23| ) being otherwise 
convergent with a logarithmic behavior. 
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4 One— loop resummation 



The success in the computation of the counterterm function K(e, a s ) is encourag- 
ing in view of a more explicit evaluation of the full form factor. It is particularly 
tempting to study the behavior of the integrals in Eq. ( |2.10| ) in the small-£ 2 region, 
where one expects to find obstructions due to the Landau pole, which should how- 
ever cancel in the ratio of the timelike to the spacelike factor, Eq. ( |2.11|) , which is 
ultraviolet dominated. 

A natural place to start is a more precise characterization of the Landau pole 
in the context of dimensional regularization. The one-loop running coupling for 
finite e can still be characterized by a mass scale, defined as the value of the renor- 
malization scale where the coupling diverges, just as in conventional dimensional 
transmutation. In fact, Eq. (|2.7|) still has a simple pole located at 



V b a s (Q 2 ; , 
and one can use the scale A 2 to define the coupling, as 



A-TTf 

<*.« 2 ) = , r , n „, n • (4-2) 



ffl - 1 



By inspection, Eqs. (|2.7|) , ( |4.1|) and ( |4.2| ) reduce to their four-dimensional counter- 
parts as e — > 0. There are however two important observations to be made. 



The running coupling defined by Eq. ( |4.2| ) vanishes as Q 2 — > for fixed e, 
provided Re e < 0, in agreement with Eq. ( |2.7| ) and with Eq. ( |2.9| ). This 
behavior for small Q 2 is simply due to dimensional counting: in fact, since we 
are keeping the regulator in place, we are studying the renormalized coupling 
9r = 9n^~ e for fixed bare coupling, in d > 4, where it vanishes as a power of 
the scale. 

For general (complex) e, the Landau pole in Eq. ( |4.1|) is located at complex 
values of the scale fi 2 . In particular, for real negative e, the location of the pole 
has a nonvanishing imaginary part provided e < — boa s (Q 2 ) / (An) , e ^ —1/n. 
Thus, in general, the Landau pole in the dimensionally regulated exponent of 
Eq. ( |2.10| ) is not on the integration contour. One expects, and I will verify 
below, that the integrals appearing in Eq. Q2.10| ) may be evaluated explicitly 
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in terms of analytic functions of a s and e. When e becomes close to 0, the 
pole migrates to the integration contour, and correspondingly the integral will 
develop a cut. Nonperturbative information (if any) will be encoded in the 
structure of the singularities at the cut, and in the corresponding higher order 
generalizations. 

Armed with these obervations, one can proceed to verify the integrability of 
logr(Q 2 ), starting with the simplest situation, in which only one-loop information 
is retained. The integration can be performed using three different methods, which 
highlight different features of the exponent. 



4.1 Integration by series 

The idea of re-expanding the one-loop running coupling, Eq. (|2.7|) , in powers of 
a s (Q 2 ) is perhaps the first that comes to mind, in view of its frequent usage in the 
context of the study of power corrections. One writes 



a 



US - 1 



(4.3) 



n=0 

and one integrates the series term by term, looking for possible factorial growth 
in the coefficients of the series obtained after integration. This would then be 
interpreted by standard methods as a sign of nonperturbative contributions. In the 
present case, changing variables in each integral in Eq. (|2.10|) as 

A 2 



A 2 -^= [£) -1, (4.4) 



one finds that all integrals as easily performed. The integration of the anomalous 
dimension 7^ generates terms that are indipendent of £, the integration variable 
of the outer integral in Eq. ( |2.10| ). These terms are divergent at the lower limit 
of integration, £ 2 = 0, but they are cancelled by the ^-independent terms in the 
expansion of the counterterm function K(e, a s ), a fact that can be explicitly verified 
to all orders by using the results of Section [| Considering for simplicity the spacelike 
form factor (obtained from Eq. ( |2.10| ) by changing Q 2 —>■ —Q 2 ), the result is then 

f-Q 2 , 9, \ 2 , 9x ~ ifa(u 2 )\ n (-l) n+1 , s 
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2e 2 ^ 




ra+1 




where I defined 

and I wrote the function G perturbatively as 

oo 

G(-l,S,e)^GW(6) 

n=l 



(4.6) 



(4.7) 



with G«(e) = C F (3 - e(C(2) - 8))/2 + C(e 2 ). 

It is apparent that the series in Eq. (|4.5| ) is far from pathological. It does not 
exhibit any factorial growth, it is alternating in sign, and definitely has a finite 
radius of convergence. In fact, one recognizes a logarithmic series in the terms with 
a single sum, whereas the terms in the nested sum are finite harmonic sums of the 
kind encountered performing Mellin transforms of Altarelli-Parisi kernels, which 
are related to polylogarithms. 

Having found that the series arising from the integration of Eq. (|2.10|) is well- 
behaved, clearly the simplest method to sum it is to perform the integrals directly, 
using Eq. Q2.7D instead of the Taylor expansion in Eq. 



4.2 Analytic expression 

Direct integration proceeds through steps that are closely related to the ones out- 
lined in the previous subsection. Through the change of variable in Eq. ( |4.4j ) one 
trivially integrates the one-loop anomalous dimension, obtaining a logarithm. The 
logarithm is not integrable over the scale £ 2 because of the singularity at £ 2 = 0, 
which is however subtracted by the contribution of the counterterm function (here 
one must use for K (e, a s ) the resummation of the leading poles, Eq. fl3.7|)). Us- 
ing once more the same change of variables, the second integration can also be 
performed. After minor reshuffling with dilogarithm identities, and inserting the 
values of 7^ and the result reads 



io g r 



a s (/i 2 



2CV 



Li; 



EL 

9\ 



a(jj, 2 ) 
a(jjL 2 ) + e 



- (7(e) In 



EL 



a(fi 2 ) 
a(fi 2 ) + e_ 
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where (7(e) = G^ (e) /Cf- Eq. (|4.8|) sums the series in Eq. ([4.5|), and has several 
rather remarkable properties. First, since the form factor is RG invariant, and I used 
the exact solution of the one-loop RG equation for a to perform the integration, 
Eq. ( |4.8|) should be independent of /i 2 to one-loop accuracy in (3{a s ). This is easily 
verified, and one obtains 



io g r 



2C F 



a s (fi 2 



i og r(-i,« s (g 2 ),e 



(4.9) 



- Li 2 



a(Q 2 



a(Q 2 



+ C(e) log 



a(Q 2 ) 



Eq. ([4.9|) , and its generalization including two-loop effects given in Section [|, are the 
central results of the present paper. The quark form factor is compactly expressed 
in terms of a simple analytic function of the coupling and of e, which is manifestly 
independent of the renormalization scale, and can be re-expanded in powers of 
®s{Q 2 ) to generate perturbative results that can be directly compared with Feynman 
diagram calculations. In particular, Eq. (|4.9| ) resums all leading (a"/e 2n ) and next- 
to-leading (a™/e 2n_1 ) infrared and collinear poles in the form factor. 

The second important feature of Eqs. ( (4.8| ) and ( |4.9j ) is that it is possible to 
study the structure of the singularities of the form factor in the complex e plane. A 
detailed analysis of this question is left for future work. Here I will simply note a few 
aspects that are apparent in Eq. (|4.9|). First, the Landau singularity, as expected, 
appears in Eq. ( |4.9| ) as one of the two branching points of a cut which is shared by 
both terms. The cut can be taken to run along the negative real axis in the e plane 
(the "physical" region), from e = —a(Q 2 ) (the Landau singularity, see Eq. ( |4.1| )) to 
e = 0. At the Landau point, both terms in Eq. (|4.9| ) diverge logarithmically. On 
the other hand, the limit e — > 0, which is the physically relevant one, can be studied 
explicitly; one finds 



logT -l,a s (Q 2 ),e 



2C f 



C(2) 



+ 



a(Q 2 



a(Q 2 



^) l0g 



'a(Q 2 



+ C(e,eloge) 



(4.10) 



Eq. ( 4.10| ) has interesting features: one observes that resumming all the leading 
poles in the logarithm of the form factor softens the singularity in a manner similar 
to what happens for the counterterm function (see Eq. (|3.7| )). In this case the 
leading singularity is a simple pole, with a residue which does not depend on the 
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coupling, an thus not on the scale either. Including higher loops in (3 corresponds to 
resumming weaker singularities, thus one might expect that it would not affect the 
leading singularity in Eq. ([4.10 ) . In Section [|, I will verify that this is indeed the case 
at two loops. Another interesting feature of Eq. ([4.1 0|) is the appearance of a term 
in the form factor of the form exp(c/a(Q 2 )), which, in the e — > limit, translates 
into a power behavior of the form (Q 2 /A 2 ) c . The term in question in the present 
case is not of direct physical significance: it leads to a positive fractional power 
of Q 2 , and it is tightly connected with the infrared divergence; in fact, as we will 



see, it cancels in the ratio (|2.11|) , which is the simplest quantity of physical interest 
that can be built out of the form factor. It is however interesting that such a term 
can arise in this way: one sees that dimensional regularization and resummation 
can combine to give well-defined (and gauge-invariant) power "corrections", which 
might be of physical interest if the present method can be generalized to construct 
infrared safe quantities to the same level of accuracy. 

As a final test of the correctness of Eq. ( |4.9| ), one can compute the ratio 

0. To the present accuracy, one finds 



T(Q 2 )/T(-Q 2 ) in the limit e 

r(l,a s (Q 2 ),e) 



log 



T(-l,a s (Q 2 ),e) 



2C f 



VK I 1 + log 



a(Q 2 



(4.11) 



a(Q 2 



- + iTr) log (l + i7ra(Q 2 ))+ 0(e) 



which agrees with Eq. (|3.7|) for the logarithmically divergent term, and agrees with 
the results of for the finite contribution, which was computed there to this same 
accuracy using Eq. (|2.11| ). Notice that the the term responsible for the power 



behavior in Q 2 has cancelled together with the simple pole. 



4.3 Integration over the coupling 



Generalizing Eqs. (4J3) and (4J3) to include two- and higher loop effects would 
appear to be nontrivial with the techniques described so far. In fact, even at two 
loops the RG equation for a cannot be solved in terms of elementary functions 
without approximations. There is, fortunately, no need to do that: one can change 
variables according to 

^ = p. , (4.12) 
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i.e. use the coupling directly as integration variable, without the need to introduce 
an explicit expression for a s (/i 2 ) in the intermediate stages. At one loop 

K = ~— \ • (4-13) 

/j, 2 a s e + f^a s 

One recognizes yet another way to characterize the Landau singularity in dimen- 
sional regularization: it arises because the one-loop f3 function with e 7^ does not 
have a (double) zero at a s = 0, as usual, but has two distinct zeroes, one at the 
origin and one located at a s = —4ire/bo. This second zero will be on the integration 
contour for real negative e, provided |e| < b a s (Q 2 )/ (47r), as expected from the ear- 
lier calculations. For e < 0, this second zero, where the (3 function is decreasing, is 
the one responsible for asymptotic freedom, while at the origin in the a s plane the 
(3 function vanishes with positive derivative, as it would in a QED-like theory, such 
as (p 4 in d = 4. In fact, this kind of behavior of dimensionally regularized theories 
(with the opposite sign of the (3 function) is familiar in statistical field theory ||18|| , 



where dimensional continuation is used, for example, to study the properties of the 
(nontrivial) 4 theory in d = 3 starting from the (trivial) case d = 4. In the present 
case, one can see from the explicit solution, Eq. (|2.7|) , that below the critical point 
e = — b a s (Q 2 ) / (4ti) the coupling decreases smoothly to 0, starting from the bound- 
ary value a s (Q 2 ). Above the critical point, on the other hand (and thus getting 
closer to the physical value e = 0), the coupling develops a Landau pole, so that it 
becomes impossible to evolve smoothly along real values of the scale. 

Inplementing the change of variables Eq. ([4.12 ), one sees that the only remaining 
dependence on the scale /j, 2 is in the upper limit of integration of the anomalous di- 
mension integral, which is now ct s (/i 2 ), and in the counterterm function K(e, a s (/x 2 )). 
These two dipendences must (and do) cancel by RG invariance, in the same way 
in which the singularities at £ 2 — > cancelled in the previous versions of the cal- 
culation. At the lower limit of integration, £ 2 = 0, one must use a s (0) = 0, in 
accordance with the arguments given above and at the beginning of Section [|; one 
then recovers directly Eq. ( |4.9| ). This third procedure of integration is clearly, and 
by far, the most straightforward; furthermore, it is readily generalizable to to more 
than one loop, as I will now show. 
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5 Two— loop resummation 



The generalization of the procedure outlined in Section [O] to include higher orders 
in the (5 function is straightforward. At I loops, /3(a s ) is a polynomial of degree Z + l, 
thus Eq. (|2.10|) will be expressed in terms of a double integral of a combination of 
rational function of a s , which in general will be computable by partial fractioning 
in terms of combinations of polylogarithms. The result will be expressed in terms 
of the perturbative coefficients of the functions fl(a s ), 7_^(a s ) and G(l, a, e), which 
are known to two loops; therefore I will briefly illustrate the procedure and the 
results in that case, where all the ingredients of the final expression are explicitly 
known, and the size remains manageable. 

At two loops one must use, for each integral over the scale of the coupling in 
Eq. (|2.10|) , the change of variables 

dfi 2 da s 1 

Partial fractioning of the denominators is achieved by introducing the two nontrivial 
zeroes of the two-loop (3 function, or, for convenience, the two zeroes of the equation 
bia 2 + boa + 4e = (with a = a s /ir), 



^--kH 1 -**)- (5 - 2) 

Note that it will be possible to verify that the result reduces to the correct limit 
at one loop by using the fact that, as b\ — > 0, a_ becomes the one-loop zero, 
a_ — > — 4e/&o, while a + diverges as a + — > —b /bi. 

The cancellation of the dependence on a s (n 2 ) between the counterterm function 
K and the anomalous dimension integral takes place as usual. The final result can 
be written in the form 



bi(a + — aJ) [ ^ v ' T v ' ) ° ^ 7ia + 



logr (-l,a s (Q 2 ),e) = ~ hJn 2 _ n S I (G {1) (e) + a + G^(e)) log (l 
+ 2 (^+aWP 



4e \ ixa + J 2bia + (a + — a_) y ira^ 
1 (f™^a j m\_ u J a + \ 



&i<2_(a + — a_) y y 7r(a + — a J) J \a + — a_ J 

+ (a + <-> a_) . (5.3) 



+ , og f 1 _^, og MQ 2 ) — 



ira + J \ ir(a + — a_) 
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It is fairly straightforward to check that Eq. ( [5.3| ) reduces to Eq. ( |4.9| ) when one 
takes {bi,^\ G^(e)} — > 0. On the other hand, the analytic structure of Eq. ( |5.3| ) 
is considerably more complicated than the one briefly discussed at one loop. In 
particular, the presence of a second nontrivial zero in (3{a s ) intoduces a new branch 
point, so that the structure of cuts becomes more complicated. One does not 
expect, however, the behavior physical quantities related to the form factor to be 
qualitatively affected by the inclusion of two-loop effects, since these are understood 
to modify the running of the coupling by an amount which is a logarithmically 
suppressed and smooth funcion of the energy. 

It is interesting to study the behavior of the two-loop resummed form factor, 
Eq. ( |5.3|) , in the neighborhood of the "physical" limit, e — > 0. One finds that, as 
expected, the simple pole in the logarithm of the form factor given in Eq. (|4.10|) 
is not affected by the inclusion of two-loop effects. The logarithmic singularities, 
however, are enhanced to log 2 strength, and also the constant terms in the limit 
e — ► receive new nontrivial contributions. One finds 



logT {-l,a s {Q ),e) = + , 3 log 

x ' o e o 



Aire ( a s {Q 2 )bi 




a s (Q 2 )b bl 1 * 



b a s (Q 2 ) V tt&o 
Aire f i + a s (Q 2 )b 1 



b a s (Q 2 ) \ nb 



o 



% \nb + ct s (Q 2 )&i/ \a s (Q 2 )b bi b x J \ nb 

+ ^wm 2C(2)) - (1 - c(2,) "f- + ° {f - e log E) • (5 - 4) 

where G$ = lim^o (e) . Also in this limit, it is easy to check that taking 
{b u ^\G^(e)} -> one gets back Eq. ( ^IQ|) . 

Although it is not proved here, one can conjecture rather confidently that the 
simple pole in \ogT(Q 2 ) is completely determined at one loop, and receives no 
corrections from higher orders in (3. The presence of the pole is in fact tightly 
connected with the one-loop zero of the f3 function, which is located at a distance 
O(e) from the origin. All higher-order zeroes have locations in the a s plane that 
are independent of e in the small e limit, thus they are not expected to introduce 
further singular behavior in this region. Logarithmic singularities, on the other 
hand, arise both from the anomalous dimension integral and from the integration 
of the function G, thus one must expect that they get higher order corrections. 
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6 Summary and perspectives 



I have discussed the quark form factor in the context of dimensionally regulated 
perturbative QCD. Building upon the results of Ref. 0, as well as upon the under- 
standing developed previously by many authors I have been able to derive 
resummed analytic expression for the form factor that are independent of the choice 
of renormalization scale, and can be systematically improved upon by including 
higher orders in the QCD (3 function. 

The main ingredient in the derivation is the consistent usage of dimensional 
regularization as analytic continuation to general complex values of the space-time 
dimension, d = 4 — 2e. It turns out that this continuation does not only regulate 
ultraviolet (as well as infrared) singularities in Feynman integrals: it also regulates 
singularities such as the Landau pole, which arise upon resummation of classes 
of Feynman diagrams, and are usually interpreted as signals of nonperturbative 
physics. How this happens is well understood, and this understanding has been 
applied as a technical tool in statistical field theory for a long time: for d > 4, 
the QCD (3 function develops an asymptotically free fixed point at a distance (9(e) 
from the origin in the a s plane, while the origin itself becomes infrared free; the 
coupling vanishes like a power of the renormalization scale, and it does so smoothly 
for sufficiently large d; when d gets close to the physical value d = 4 the coupling 
develops a Landau pole for real values of the scale; this pole is then found on 
the integration contour of resummed expressions for QCD amplitudes, which, as a 
consequence, develop cuts with computable branching points and discontinuities. 

The versatility of dimensional regularization was exploited here also to derive 
compact resummed analytic expressions for the counterterm function K(e, a s ): the 
fact that the renormalization group equation for a pure counterterm (with a finite 
anomalous dimension) turns into a recursion relation is well known; the fact that for 
the quark form factor the recursion relation is completely solvable to all orders in 
/3, and the resulting towers of poles can be summed, provided a strong motivation 
to pursue the calculation for the complete form factor. 

The calculations presented in this paper are not directly applicable to physical 
processes: the paper, after all, deals with a divergent quantity^. There are however 
several lines of enquiry that are opened by these results, and I believe are worth 



3 The modulus squared of Eq. ( 4.1l| ) is an exception, as it enters directly the resummed Drell- 
Yan cross section; it was, however, computed already in Ref. (]]] 
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pursuing. 

First of all, it would be of great interest to extend results such as Eq. ( |4.9| ) 
and Eq. (|5.3| ) to more general QCD amplitudes, in particular to amplitudes with 
more than two colored particles. This is not straightforward, because more com- 
plicated amplitudes have a nontrivial tensor structure in color space, and evolution 
equations like Eq. ( |2.2|) correspondingly turn into matrix equations. It should be 
emphasized however that much of the necessary technology has already been devel- 
oped in the context of the resummation of large-x logarithms for the production 
of heavy quark pairs and other colored final states in hadronic collisions [fTO] . In 
that case one resums large corrections due to the presence of + distributions at 
each order in perturbation theory, arising from the imperfect cancellation of real 
and virtual contributions. In the present case one would be interested in purely 
virtual corrections contributing to the cross section with terms proportional to 5 
functions. I believe that a generalization of the present results in this direction 
should be possible, and work to achieve it is under way. 

A more ambitious goal would be to attempt the construction of finite, physical 
quantities to the same accuracy achieved here for the form factor. This would re- 
quire controlling the infrared/collinear singularities of dimensionally regulated real 
emission amplitudes to all orders. While this is not likely to be possible in a general 
multi-scale process, progress might conceivably be made in the simplest cases, such 
as the inclusive Drell-Yan cross section, by introducing suitable approximations to 
the <i-dimensional cross section. 

On a more speculative note, one may wonder whether the analytic structure 
uncovered in the present paper may shed some light on nonperturbative features 
of QCD amplitudes. It was already noticed in the comments to Eq. (|4.10 ) that in 
the neighborhood of e = one finds terms that translate into powers of the ratio 
Q 2 /A 2 . More generally, the analytic features of the form factor near the Landau 
cut might provide suggestions to build models for the nonperturbative behavior of 
simple QCD amplitudes, as well as for the behavior of the coupling itself, much as 
renormalon calculations did |T2|, |H| . 

In the context of finite order perturbative QCD calculations, the existence of 
genuinely i?G-invariant expressions such as Eqs. ( |4.9[ ) and Q5.3p might prove useful 
as a model for the renormalization scale dependence of more complicated QCD 
amplitudes. Not surprisingly, the results for the form factor are naturally expressed 
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in terms of a s (Q 2 ), since Q 2 is the only perturbative scale in the amplitude. One 
can however go back to equations such as ( |4.8| ), choose a different renormalization 
scale, say \i = Q/2, re-expand the form factor to a finite perturbative order, and 
study the variations of the answer with \i. 

Finally, it should be mentioned that large corrections to perturbative QCD 
amplitudes due to analytic continuation (the "7r 2 terms" ) do not arise only in quark- 
initiated process. A case in point is Higgs production via gluon fusion, in an effective 



theory in which the top quark has been integrated out |[20|| . Also in that case the 
present techniques may well prove useful and provide a very accurate and rigorous 
evaluation of the impact of these terms to all orders in the coupling. 

In summary, I believe that the results presented here should be of interest for 
the study of several aspects, both formal and practical, of perturbative QCD. Much 
work remains to be done. 
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